skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Yarza, Ricardo"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract Once per ≈104–105yr, an unlucky star may experience a close encounter with a supermassive black hole (SMBH), partially or fully tearing apart the star in an exceedingly brief, bright interaction called a tidal disruption event (TDE). Remnants of partial TDEs are expected to be plentiful in our Galactic center, where at least six unexplained, diffuse, star-like “G objects” have already been detected, which may have formed via interactions between stars and the SMBH. Using numerical simulations, this work aims to identify the characteristics of TDE remnants. We take 3D hydrodynamic FLASH models of partially disrupted stars and map them into the 1D stellar evolution code MESA to examine the properties of these remnants from tens to billions of years after the TDE. The remnants initially exhibit a brief, highly luminous phase, followed by an extended cooling period as they return to stable hydrogen burning. During the initial stage (≲105yr) their luminosities increase by orders of magnitude, making them intriguing candidates to explain a fraction of the mysterious G objects. Notably, mild TDEs are the most common, and result in the brightest remnants during this initial phase. However, most remnants exist in a long-lived stage where they are only modestly offset in temperature and luminosity compared to main-sequence stars of equivalent mass. Nonetheless, our results indicate remnants will sustain abnormal, metal-enriched envelopes that may be discernible through spectroscopic analysis. Identifying TDE survivors within the Milky Way could further illuminate some of the most gravitationally intense encounters in the Universe. 
    more » « less
    Free, publicly-accessible full text available August 27, 2026
  2. Abstract Stars grazing supermassive black holes on bound orbits may produce periodic flares over many passages, known as repeating partial tidal disruption events (TDEs). Here, we present 3D hydrodynamic simulations of Sun-like stars over multiple tidal encounters. The star is significantly restructured and becomes less concentrated as a result of mass loss and tidal heating. The vulnerability to mass loss depends sensitively on the stellar density structure, and the strong correlation between the fractional mass loss ΔM/M*and the ratio of the central and average density ρ c / ρ ¯ , which was initially derived in disruption simulations of main-sequence stars, also applies for stars strongly reshaped by tides. Over multiple orbits, the star loses progressively more mass in each encounter and is doomed to a complete disruption. Throughout its lifetime, the star may produce numerous weak flares (depending on the initial impact parameter), followed by a couple of luminous flares whose brightness increases exponentially. Flux-limited surveys are heavily biased toward the brightest flares, which may appear similar to the flare produced by the same star undergoing a full disruption on its first tidal encounter. This places new challenges on constraining the intrinsic TDE rates, which need to take repeating TDEs into account. Other types of stars with different initial density structures (e.g., evolved stars with massive cores) follow distinct evolution tracks, which might explain the diversity of the long-term luminosity evolution seen in recently uncovered repeaters. 
    more » « less
    Free, publicly-accessible full text available January 16, 2026
  3. Abstract We present a start-to-end simulation aimed at studying the long-term fate of high-mass X-ray binaries and whether a Thorne–Żytkow object (TŻO) might ultimately be assembled. We analyze results from a 3D hydrodynamical simulation that models the eventual fate of LMC X-4, a compact high-mass X-ray binary system, after the primary fills its Roche lobe and engulfs the neutron star companion. We discuss the outcome of this engulfment within the standard paradigm of TŻO formation. The post-merger angular momentum content of the stellar core is a key ingredient, as even a small amount of rotation can break spherical symmetry and produce a centrifugally supported accretion disk. Our findings suggest the inspiraling neutron star, upon merging with the core, can accrete efficiently via a disk at high rates (≈10−2Ms−1), subsequently collapsing into a black hole and triggering a bright transient with a luminosity and duration typical of an ultra-long gamma-ray burst. We propose that the canonical framework for TŻO formation via common envelope needs to be revised, as the significant post-merger accretion feedback will unavoidably unbind the vast majority of the surrounding envelope. 
    more » « less
  4. Abstract Triple stellar systems allow us to study stellar processes that cannot be attained in binary stars. The evolutionary phases in which the stellar members undergo mass exchanges can alter the hierarchical layout of these systems. Yet, the lack of a self-consistent treatment of common-envelope (CE) in triple-star systems hinders the comprehensive understanding of their long-term fate. This paper examines the conditions predicted around binaries embedded within CEs using local 3D hydrodynamical simulations. We explore varying the initial binary separation, the flow Mach number, and the background stellar density gradients as informed by a wide array of CE conditions, including those invoked to explain the formation of the triple system hosting PSR J0337+1715. We find that the stellar density gradient governs the gaseous drag force, which determines the final configuration of the embedded binary. We observe a comparable net drag force on the center of mass but an overall reduction in the accretion rate of the binary compared to the single-object case. We find that, for most CE conditions, and in contrast to the uniform background density case, the binary orbital separation increases with time, softening the binary and preventing it from subsequently merging. We conclude that binaries spiraling within CEs become more vulnerable to disruption by tidal interactions. This can have profound implications on the final outcomes of triple-star systems. 
    more » « less
  5. Common-envelope evolution is a stage in binary system evolution in which a giant star engulfs a companion. The standard energy formalism is an analytical framework to estimate the amount of energy transferred from the companion's shrinking orbit into the envelope of the star that engulfed it. We show analytically that this energy transfer is larger than predicted by the standard formalism. As the orbit of the companion shrinks, the mass it encloses becomes smaller, and the companion is less bound than if the enclosed mass had remained constant. Therefore, more energy must be transferred to the envelope for the orbit to shrink further. We derive a revised energy formalism that accounts for this effect, and discuss its consequences in two contexts: the formation of neutron star binaries, and the engulfment of planets and brown dwarfs by their host stars. The companion mass required to eject the stellar envelope is smaller by up to 50% , leading to differences in common-envelope evolution outcomes. The energy deposition in the outer envelope of the star, which is related to the transient luminosity and duration, is up to a factor of ≈7 higher. Common-envelope efficiency values above unity, as defined in the literature, are thus not necessarily unphysical, and result at least partly from an incomplete description of the energy deposition. The revised energy formalism presented here can improve our understanding of stellar merger and common-envelope observations and simulations. 
    more » « less
  6. Abstract Stars grazing supermassive black holes (SMBHs) on bound orbits may survive tidal disruption, causing periodic flares. Inspired by the recent discovery of the periodic nuclear transient ASASSN-14ko, a promising candidate for a repeating tidal disruption event (TDE), we study the tidal deformation of stars approaching SMBHs on eccentric orbits. With both analytical and hydrodynamic methods, we show the overall tidal deformation of a star is similar to that in a parabolic orbit provided that the eccentricity is above a critical value. This allows one to make use of existing simulation libraries from parabolic encounters to calculate the mass fallback rate in eccentric TDEs. We find the flare structures of eccentric TDEs show a complicated dependence on both the SMBH mass and the orbital period. For stars orbiting SMBHs with relatively short periods, we predict significantly shorter-lived duration flares than those in parabolic TDEs, which can be used to predict repeating events if the mass of the SMBH can be independently measured. Using an adiabatic mass-loss model, we study the flare evolution over multiple passages, and show the evolved stars can survive many more passages than main-sequence stars. We apply this theoretical framework to the repeating TDE candidate ASASSN-14ko and suggest that its recurrent flares originate from a moderately massive ( M ≳ 1 M ⊙ ), extended (likely ≈10 R ⊙ ), evolved star on a grazing, bound orbit around the SMBH. Future hydrodynamic simulations of multiple tidal interactions will enable realistic models on the individual flare structure and the evolution over multiple flares. 
    more » « less
  7. Abstract The nearby type II supernova, SN 2023ixf in M101 exhibits signatures of early time interaction with circumstellar material in the first week postexplosion. This material may be the consequence of prior mass loss suffered by the progenitor, which possibly manifested in the form of a detectable presupernova outburst. We present an analysis of long-baseline preexplosion photometric data in theg,w,r,i,z, andyfilters from Pan-STARRS as part of the Young Supernova Experiment, spanning ∼5000 days. We find no significant detections in the Pan-STARRS preexplosion light curves. We train a multilayer perceptron neural network to classify presupernova outbursts. We find no evidence of eruptive presupernova activity to a limiting absolute magnitude of −7 mag. The limiting magnitudes from the full set ofgwrizy(average absolute magnitude ≈ −8 mag) data are consistent with previous preexplosion studies. We use deep photometry from the literature to constrain the progenitor of SN 2023ixf, finding that these data are consistent with a dusty red supergiant progenitor with luminosity log L / L ≈ 5.12 and temperature ≈ 3950 K, corresponding to a mass of 14–20M
    more » « less
  8. The engulfment of substellar bodies (SBs) such as brown dwarfs and planets has been invoked as a possible explanation for the presence of SBs orbiting subdwarfs and white dwarfs, rapidly rotating giants, and lithium-rich giants. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an SB engulfed in a stellar envelope. We model the SB as a rigid body with a reflective boundary because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory during engulfment. We find that SB engulfment can increase the stellar luminosity of a 1M⊙ star by up to a few orders of magnitude for timescales of up to a few thousand years when the star is ≈10R⊙ and up to a few decades at the tip of the red giant branch. We find that no SBs can eject the envelope of a 1M⊙ star before it evolves to ≈10R⊙ . In contrast, SBs as small as ≈10MJup can eject the envelope at the tip of the red giant branch, shrinking their orbits by several orders of magnitude in the process. The numerical framework we introduce here can be used to study the dynamics of planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB. 
    more » « less
  9. The engulfment of substellar bodies (SBs) such as brown dwarfs and planets has been invoked as a possible explanation for the presence of SBs orbiting subdwarfs and white dwarfs, rapidly rotating giants, and lithium-rich giants. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an SB engulfed in a stellar envelope. We model the SB as a rigid body with a reflective boundary because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory during engulfment. We find that SB engulfment can increase the stellar luminosity of a 1M⊙ star by up to a few orders of magnitude for timescales of up to a few thousand years when the star is ≈10R⊙ and up to a few decades at the tip of the red giant branch. We find that no SBs can eject the envelope of a 1M⊙ star before it evolves to ≈10R⊙ . In contrast, SBs as small as ≈10MJup can eject the envelope at the tip of the red giant branch, shrinking their orbits by several orders of magnitude in the process. The numerical framework we introduce here can be used to study the dynamics of planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB. 
    more » « less
  10. Abstract Young terrestrial worlds are critical test beds to constrain prevailing theories of planetary formation and evolution. We present the discovery of HD 63433 d—a nearby (22 pc), Earth-sized planet transiting a young Sun-like star (TOI-1726, HD 63433). HD 63433 d is the third planet detected in this multiplanet system. The kinematic, rotational, and abundance properties of the host star indicate that it belongs to the young (414 ± 23 Myr) Ursa Major moving group, whose membership we update using new data from the third data release of the Gaia mission and TESS. Our transit analysis of the TESS light curves indicates that HD 63433 d has a radius of 1.1Rand closely orbits its host star with a period of 4.2 days. To date, HD 63433 d is the smallest confirmed exoplanet with an age less than 500 Myr, and the nearest young Earth-sized planet. Furthermore, the apparent brightness of the stellar host (V≃ 6.9 mag) makes this transiting multiplanet system favorable to further investigations, including spectroscopic follow-up to probe the atmospheric loss in a young Earth-sized world. 
    more » « less